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<^ ; Abstract 



. The observation of scaling in processes in which a weakly interacting probe 

. delivers large momentum q to a many-body system reflects the dominance of 

incoherent scattering off target constituents. While a suitably defined scaling 
function can provide rich information on the internal dynamics of the target, 
in general its extraction from the measured cross section requires careful con- 

"O . 

sideration of the nature of the interaction driving the scattering process. The 



analysis of deep inelastic electron-proton scattering in the target rest frame 
within standard many-body theory naturally leads to the emergence of a scal- 
ing function that, unlike the commonly used structure functions F± and Fz, 
can be identified with the intrinsic proton response. The implications for the 
theoretical analysis of deep inelastic scattering and the interpretation of the 
data are discussed. 

PACS numbers: 13.60.Hb,24.10.Cn,61.12.Bt 



1 



Typeset using REVTgX 



I. INTRODUCTION 



Scaling is observed in a variety of scattering processes involving many-body systems 
[IJ. For example, at large momentum transfer |q| the response of liquid helium measured 
by inclusive scattering of thermal neutrons, which in general depends upon both q and the 
energy transfer u, exhibits a striking scaling behavior, i.e. it becomes a function of the single 
variable y = (m/\q\)(u — q 2 /2m), m being the mass of the helium atom M. Scaling in a 
similar variable occurs in inclusive electron-nucleus scattering at |q| > 500 MeV and electron 
energy loss v < Q 2 /2M, where Q 2 = |q| 2 — v 2 and M is the nucleon mass Another most 
celebrated example is scaling of the deep inelastic proton structure functions, measured by 
lepton scattering at large Q 2 , in the Bjorken variable x = Q 2 /2Mu [|J. 

The relation between Bjorken scaling and scaling in the variable y, whose definition and 
interpretation emerge in a most natural fashion from the treatment of the scattering process 
in the target rest frame within many-body theory, has been discussed by many authors (see, 
e.g., ref. |J). Recently, deep inelastic data have been also shown to scale in the variable 
y — v — |q| ||, related to both y and the Nachtman variable £ J7|, which in turn coincides 
with x in the Q 2 —>■ oo limit. 

The observation of scaling in processes driven by different interactions clearly indicates 
that its occurrence reflects the dominance of a common reaction mechanism, independent of 
the underlying dynamics. In all instances, scaling is indeed a consequence of the onset of the 
kinematical regime in which scattering of a weakly interacting probe by a composite target 
reduces to the incoherent sum of elementary scattering processes involving its constituents. 

While the primary goal of scaling analysis is the identification of the dominant reaction 
mechanism, it has to be also emphasized that the scaling variable has a straightforward 
physical interpretation, and a suitably defined scaling function, being directly related to the 
target response, contains a great deal of dynamical information. 

In general, extracting the target response from the mesasured cross section requires care- 
ful consideration of the nature of the interaction driving the scattering process. In neutron- 
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liquid helium scattering, due to the scalar character of the probe-constituent coupling, the 
cross section coincides with the response up to a kinematical factor 0. On the other hand, 
in electromagnetic processes the internal dymanics of the target and the probe-constituent 
interaction are nontrivially coupled. As a consequence, to obtain the target response in the 
case of electron-nucleus scattering one has to devise a procedure to remove from the data 
the dependence upon the elementary electron-nucleon cross section ||. In this paper I will 
discuss the extension of this method to the analysis of deep inelastic scattering (DIS) of 
electrons by protons. 

The theoretical treatment of scattering off a many-body system and the emergence of 
scaling as a consequence of the assumptions involved in the impulse approximation (IA) are 
outlined in Section II, where the case of neutron-liquid helium scattering is considered as a 
pedagogical example. The procedure leading to identify the target response and the scaling 
variable in the more complex case of electron-nucleus scattering is traced in Section III, while 
Section IV describes the application of the same procedure to deep inelactic electron-proton 
scattering, and compares the resulting analysis to the standard Bjorken scaling analysis. 
The implications of the proposed approach for both the theoretical description of DIS and 
the interpretation of the data are discussed in Section V. Finally, Section VI summarizes 
the main results and states the conclusions. 

II. SCATTERING OFF MANY-BODY SYSTEMS IN THE IA REGIME AND 

y-SCALING 

A. Target response function 

Let us consider scattering off a norelativistic bound system consisting of N pointilke 
particles of mass m, and assume that the probe-target interaction be scalar and weak, so 
that Born approximation can safely be used. The differential cross section of the process 
in which a beam particle, carrying momentum k and energy E, is scattered into the solid 
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angle dfl with energy E' = E — v and momentum k' can be written 

do a Ik' 



d£ldE' 4vr Ikl 



s(a,v), (i) 



where a is the probe-constituent total cross section. The response function 5(q, u), contain- 
ing all the information on the structure of the target, is defined as 

S(q, i/) = £ I (n|ft,|0> \ 2 5(u + E - E n ) (2) 

n 

r dt 

= /^e^ (0|pt(t) Pq ( )|0> , (3) 

where 

p q (t) = e iHt Pq ^ im , (4) 
if is the hamiltonian describing the internal dynamics of the target and 

Pq = Yj a P+q fl P > ( 5 ) 

p 

aj, and a p being constituent creation and annihilation operators, respectively. The target 
ground and final states |0) and \n) satisfy the many-body Schrodinger equations H\0) = 
Eq\0) and H\n) = E n \n). 

Rewriting Eq.(0) in coordinate space we obtain 



5(q, u) = E 



dR {n\R)Y,^' Tl W) 



S(u + E + E n ), (6) 



i=i 

where R = (r±, . . . , i*n) specifies the target configuration, (i?|0) and (i?|n) denote the tar- 
get initial and final state wave functions, respectively, and the index i labels the struck 
constituent. 

The main assumption underlying IA is that, as the space resolution of a probe delivering 
momentum q is ~ l/|q|, at large enough |q| (typically |q| 3> 2ir/d, d being the average 
separation between target constituents) the target is seen by the probe as a collection of 
individual particles. In addition, final state interactions (FSI) between the hit constituent, 
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carrying large momentum ~ q, and the residual (N — l)-particle system are assumed to be 
negligibly small. 

In the IA regime the scattering process reduces to the incoherent sum of elementary 
processes involving only one constituent, the remainig (N — 1) particles acting as spectators, 
and Eq.(0) takes the simple form 

2 5(u + E + E n ). (7) 

On account of the absence of FSI, the IA final state |n), carrying total momentum q, exhibits 
the product structure 

|n) = |p , ,^) = |p')®|^(q-p')) , (8) 

its energy being given by 

E n = E p , + E n . (9) 

In the above equations |p') denotes the state describing a free particle carrying momentum 
p' and energy E p > = |p'| 2 /2m, while En is the energy of the the spectator (N — l)-particle 
system, left in the state \JZ) with momentum q — p'. As a consequence, the sum over final 
states appearing in Eq.(f|) can be carried out replacing 

J»(n| - / dV |p'>(p'| £ |^(q-p'))(^(q-p')l • (10) 

Substitution of Eqs.(|3D-(|T0D into Eq.(^) leads to (see, e.g., ref. [|J) 

5(q, u) = J d 3 p J de P(p, e)8(u + e - E p+q ) , (11) 

where the function 

P(p, e) = Y / \(0\P,K)\ 2 5(e + E n - E ) (12) 
n 

gives the probability of finding a target constituent with momentum p and energy e in the 
target ground state. 
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%^) = EE JdR (n\R) 



?qr, 



(J2|0> 



B. y-scaling: neutron scattering off liquid helium as an example 



Eqs.flTTD and (|T2"D show that the IA response only depends upon q and v through the 
energy-conserving ^-function, requiring 

v + E - E n - £ p+q = . (13) 

The occurrence of scaling, i.e. the fact that, up to a kinematical factor K(\q\,u), S*(q, v) 
becomes a function of a single variable, simply reflects the fact that in the IA regime, in which 
energy conservation is expressed by Eq.(|l3|), q and v are no longer independent variables. 
One can then define a new variable y = y(q, v) such that, as |q| — > oo, 

K(\tilv)S(^v)^F(y) . (14) 

All the above results can be immediatly applied to the case of neutron scattering off 
liquid helim, in which the energy dependence of P(p, e) can be safely neglected and Eq.(|13D 
takes the form || 



m being the mass of the helium atom. It follows that, definin 

|2\ 

|q| y 2m 
one can rewrite the response in the form 



S(q,v) = f- 27r \p\d\p\ n(|p|) , (17) 

|q| J\v\ 

with the constituent momentum distribution defined as 

n(|p|) = JdeP{p,e). (18) 
Eq.(|l7D implies that in the |q| — >• oo limit 

M S(q, v) - F(y) , (19) 



the kinematical factor being given by i^(|q|, v) = \du/dpu\ = \q\/m, with p\\ = p ■ q/| q| . 

Fig.0 displays the behavior of F(y), defined by the above equation, measured in neutron 
scattering off superfluid 4 He at T =1.6 °K 0. It clearly appears that the curves corre- 
sponding to |q| > 15 A -1 lie on top of one another, indicating the onset of the scaling 
regime. 

The variable y defined by Eq. fll6|) does have a straightforward physical interpretation, as 
does the scaling function of Eq . (|19D . The scaling variable can be identified with the initial 
longitudinal momentum of the struck atom, p», while F(y) can be related to the momentum 
distribution through 

n ,n 11 (dF\ 

The relation linking the response in the scaling regime to w(|p|) has been been extensively 
exploited to extract momentum distributions of normal and superfluid 4 He from neutron 



scattering data [10 . 



III. SCALING IN ELECTRON SCATTERING 

A. Scattering cross section in the IA regime 

The unpolarized electron scattering cross section is usually written in the form (see, e.g., 
ref. H) 

B L„W , (21) 



dVLdE' Q 4 E 

where E and E' denote the initial and final electron energy, respectively, Q 2 = —q 2 , q = 
k — k' = (v, q) being the four-momentum transfer, and a is the fine structure constant. The 
electron tensor L^ u is fully specified by the measured kinematical variables. In the limit of 
ultrarelativistic electrons it reads 



— 2 



kX + Kk'-g^(kk') . (22) 
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The information on the structure of the target is contained in the tensor W^, which can be 
written in a form reminiscent of Eqs.(Q) and ([!]) replacing the density fluctuation operator 
p q with the electromagnetic current J 7 '. The resulting expression reads 

= ^(0|^V)^l^|0)^(Po + q ~ P n ) (23) 

n 

= J-0^e tqx (0\r\x)r(0)\0), (24) 

where P = (M, 0), M being the target mass, and P n denote the four-momenta of the target 
initial and final states, respectively. 

In the IA scheme the target tensor of Eqs. fl2"3] ) and fl2"3| ) is replaced by a weighted sum of 
tensors describing the electromagnetic structure of target constituents (compare to Eq. (|TTT) ) : 

W ^^J2 J d4 P P (P) q)6(v + p - £ p+q ) , (25) 

i 

where E p+Cl = \J\p + q| 2 + m 2 , m being the constituent mass, p = (E p ,p), with E p = 
\J\p\ 2 + m 2 , and q = (v, q), with v = E p+Cl — E p = v + p — E p . The distribution function 
P(p) yields the probability to find a constituent carrying four-momentum p = (po,p) in 
the target ground state. It has to be emphasized that the tensor describes the electro- 
magnetic interaction of the z-th target constituent in free space. Hence, the above equation 
shows that the IA formalism allows one to describe scattering off a bound constituent in 
terms of the electromagnetic tensor associated with a free constituent, binding effects being 
taken care of by replacing p, q with p, q (see Appendix A). 

Note that, on account of the replacement q — >• q, wf is manifestly non gauge- invariant, 
as q^w^ 7^ 0. A somewhat ad hoc prescription to restore gauge invariance, widely used in 
the theoretical analysis of electron-nucleus scattering data, has been proposed by De Forest 
back in the 80s [|T^] (see Appendix B). More recently, the problem of reconciling IA and 
current conservation in DIS has been discussed in ref. ||13|| . It has to be emphasized that, 
although in general violation of gauge invariance is an unpleasant feature inherent in the 
IA scheme, it does not play a role in the |q| — ► oo limit, in which the non gauge-invariant 
contributions to the current can be shown to be vanishingly small. 
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Substitution of the IA target tensor fl25|) into Eq.([2~TD yields 
da 



dQdE 



(26) 



where (dai/dQdE') is the elementary electron-constituent cross section, whose structure will 
be discussed in the next Section. 

B. Extracting the target response from the cross section 



From Eqs.(P), (gD and (§D it follows that 



( 



da, \ _ a 2 E' 
\dttdE>) ~ Q* ~E L ^" - 



— — L u „wr = I — 



da 
dtt 



M 



a i2 + 2 a h tan - 



(27) 



where (da/dVL)M denotes the Mott cross section, 6 is the electron scattering angle and 

is the electromagnetic tensor associated with the z-th constituent, corrected to restore gauge 

invariance, whose definition is given in Appendix B. The functions a^ and <Xj 2 can be written 

,2 , 1 p 2 ± , 



m 

E p E p+<i 



Wi i + 2m> Wi > 



(28) 



and 



m 



E p E p+n |q|* \ Wh [f~ 1 ) + ^ 



Q' 



2\ E P~ U TV 



2 2 
I PI 



(29) 



iqr v ' <r j 2 

where p± is the component of the constituent momentum perpendicular to the momentum 
transfer and and Wi 2 are the constituent structure functions. 

Substitution of Eqs. (p7|) -(p9|) into Eq . (f26|) finally leads to the familiar expression of the 
electron scattering cross section in terms of the two target structure functions W% and W<i- 



d 2 a 
dVLdE' 



da 

In 



M 



W 2 + 2Wi tan 2 



e 



(30) 



with 



Wx,2 = E Jd 4 pP(p)a ll , 2 (p,q)6 



v + Po - v |p + q| + m ' 



(31) 
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Eq.(|26|) shows that in the case of electron scattering the target structure function does 
not appear as a multiplicative factor in the cross section. In principle, the p-dependence 
of the elementary cross section, coupling the internal dynamics of the target to the probe- 
constituent interaction, prevents one from extracting the structure function from the data. 
In addition, as the integrand in the right hand side of Eqs.fl2"E|) depends upon q through 
both (dai/dQdE') and the argument of the energy conserving ^-function, the occurrence of 
y-scaling can no longer be established using the simple argument of Section II. 

The above problems can be circumvented exploiting the weak momentum and energy 
dependence of the elementary cross section. Assuming that the distribution function P(p) 
be a rapidly decreasing function, the elementary cross section can be evaluated at a constant 
p = p, corresponding to the peak of P{p), and moved out of the integral. As a result, one 
can readily identify the target response with the ratio 



Note that, as S(q, v) defined by the above equation depends upon q only through the energy 
conserving delta function, it is also expected to exhibit ^/-scaling. 



The procedure described in the previous Section has been extensively employed to analyze 
electron- nucleus {eA) scattering data. Although in this case the target constituents (Z 
protons and N = A — Z neutrons) are not structureless, ^-scaling still occurs as long as the 
internal degrees of freedom of the constituents are not excited by the interaction with the 
electron probe, i.e. as long as the elementary electron-nucleon scattering process is elastic. 
In this case, the nucleon structure functions entering Eqs.fl2~8|) and (|29|) can be written 
(i = p, n) 




(32) 



C. y-scaling in electron- nucleus scattering 




(F h + nF i2 f 



(33) 



and 
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^2 



n 4m 21 1 12 j 



(34) 



where and F i2 are the Dirac and Pauli form factors, respectively, and Ki denotes the 
nucleon anomalous magnetic moment. 

In analogy with the case of neutron scattering off liquid helium, the requirement of energy 
conservation provides the definition of the scaling variable through 



v + M - sj(y + |q|) 2 + m 2 — \Jy 2 + (M — m + A)) 2 = 0, (35) 

where B is the (positive) minimum nucleon binding energy. The corresponding scaling 
function is 

F(v) = K {q ,„) (^) (36) 



where K(q,y) = \du/dp\\\ } 



p=p- 



Fig.S shows the cross sections recently measured at the Thomas Jefferson National Accel- 



erator Facility scattering 4 GeV electrons off an Iron target 14 ] . The corresponding scaling 



functions, obtained from the definitions of Eqs.(^) and (|36"D , are shown in fig.^J It clearly 



appears that the data, covering the range of momentum transfer 1 < |q| < 4 GeV, exhibit 
a striking scaling behaviour at y < 0. The scaling violations observed in the region of posi- 
tive y, corresponding to large electron energy loss, have to be ascribed to the occurrence of 
inelastic electron-nucleon scattering. 

Although the scaling variable defined by Eq.([$5|) can still be related to the longitudinal 
momentum of the struck nucleon, in the case of eA scattering the constituent binding energy 
plays a significant role, making it more difficult to establish a direct relation between F(y) 
and the nucleon momentum distribuution n(|p|) . A procedure to extract n(|p|) from eA 
data in the ^/-scaling regime taking into account binding corrections has been developed in 
ref. p|. 



IV. SCALING IN DEEP INELASTIC SCATTERING AND PROTON RESPONSE 
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A. The proton as a many-body system 



Let us now make the rather strong assumption that the proton can be viewed as a 
many-body system consisting of bound pointlike Dirac particles of mass m and charge 
(in units of the magnitude of the electron charge), and extend the analysis described in the 
previous Sections to deep inelastic electron-proton (ep) scattering. In this scenario, as in the 



standard parton model of DIS |lT|, one assumes that over the short spacetime scale relevant 
to the scattering process confinement does not play a role, so that proton constituents can 
be described in terms of physical states. It has to be emphasized, however, that, unlike the 
parton model, the present approach does not involve the additional assumption that proton 
constituents be on mass shell. 

The structure functions appropriate for the case of pointlike constituents are (compare 
to Eqs.(||) and flU) 

w - = ~ e * h (37) 



and 



(38) 



The requirement of energy conservation, expressed by the equation 

1 



v + Po - V |p + q| 2 + m 2 = v + p - p\\ - |q| + O \ — J = , (39) 



where p = M — En and En = y |p| 2 + M n , M n being the mass of the spectator system, 
implies that, as |q| — > oo, the quantity 

y = v- |q| =p|| -p (40) 

becomes independent of q. Hence, in this limit S^q, z/), defined as in eq. (|32|) , is expected to 
exhibit scaling in the variable y 0, i.e. 

S(q,i/)-F(y). (41) 
12 



Note that in this case the kinematical factor entering the definition of the scaling function 
is \dv/dp\\ | = 1. 

It has to be pointed out that y does not have the same physical interpretation as the 
variable y discussed in the previous sections: it does not coincide with the constituent 
longitudinal momentum. However, as po is independent of q, scaling in y = p» necessarily 
implies scaling in y, and viceversa. The motivation for choosing y as scaling variable in DIS 
will be discussed in the next Section. 



B. y-scaling analysys of DIS data 

According to the IA picture, the ^-scaling function can be obtained dividing either struc- 
ture function by the appropriate contribution to the elementary cross section. In fact, from 
Eqs.(gg)-(|32|) and ([41]) it follows that, in the |q| ->• oo limit, 

5(q, y) = ^ = ^ - F(y) , (42) 

where 

°i,2 = ( f u, 2 )p=p > ( 43 ) 

i 

with p = (po;Pmm)j Pmin being the minimum constituent momentum allowed in the kine- 
matics specified by q and v. The magnitude of p m i n is given by 



1 

| Pmin | 7> 



Ml-{y + Mf 



y + M 

while po, related to the mass of the spectator system through 



(44) 



Po = M-y/\p min \ 2 + M^, (45) 
can be simply parametrized in terms of the positive quantity B according to 

M n = M -m + B . (46) 
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The above discussion shows that the scaling analysis in terms of y involves two parame- 
ters: the constituent mass, m, and B , playing the role of the constituent binding energy of 
many-body theory. 

Fig.| shows the quantities Wi/Wi (upper panel) and Wijai (lower panel), obtained from 
DIS data taken at SLAC [TfJ and CERN fT7| , |T8| and rearranged in bins of constant |q| 
centered at 11, 19 and 27 GeV, plotted as a function of y. The structure functions W\ have 
been obtained from the tabulated F 2 = vWi using the parametrization of R = <jl/ 1 o~t of 
ref. [|l]|], while the elementary cross sections have been evaluated from Eqs.(|28|)-(^9|) and 
(j37|)- (|38l) , with m = 300 MeV and Bq = 200 MeV. It clearly appears that in both cases 
scaling sets in at |q| > 10 GeV. The second feature predicted by the IA analysis, i.e. that 
Wi/ai and W2/W2 scale to the same function F(y), is illustrated in fig.[5|. 

It has to be emphasized that the occurrence of y scaling does not depend upon the values 
of either m or Bq. In particular, it does not require that the constituent mass be negligibly 
small, nor that the struck constituent be on mass shell. Varying m and B in the range 
10 — 300 MeV does not appreciably affect the scaling behavior displayed in figsfj and [5|. 
However, it does affect the scaling function extracted from the data. 

While increasing m at fixed B leads to a shift of F(y) towards lower y for y > — .5 GeV, 
increasing Bq at fixed constituent mass produces a shift in the opposite direction. As a 
result, the scaling function turns out to be only sensitive to the difference m — B , i.e. to the 
mass of the recoiling spectator system M-ji- This feature is illustrated in fig.Fig.|6|, showing 
that the scaling functions corresponding torn = 110 MeV, B = 10 MeV and m = 300 MeV, 
B = 200 MeV, yielding the same value of M-n, lie on top of one another. 



C. Comparison between y- and x-scaling analysis 

DIS data are usually analyzed in terms of the two dimensionless structure functions 
F\ = MWi and F2 = vW2- In the Bjorken limit Q 2 , v — > 00, with Q 2 / v finite and v/\q\ — > 1, 
both Fi and F 2 exhibit scaling, i.e. they become functions of a single variable x = Q 2 /2Mv. 
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The first issue to be addressed to establish a connection between x- and ^-scaling is the 
relation between the scaling variables. As pointed out in ref. |J, y is trivially related to 
another variable commonly used in the context of DIS, the Nachtman variable £ 0, through 

= ^ (47 ) 

M 1 + Jl + 4M 2 x 2 /Q 2 ' 



where x is the Bjorken variable. From the above equation it follows that in the Q 2 — > oo 
limit y coincides with x up to a constant factor. Comparison between x and —y/M at 
constant |q| = 20 MeV shows that the difference is in fact less then 2 % over the whole 
< x < 1 range. 

The correspondence between y and £ makes it clear why y is better suited than y for the 
analysis of DIS. In addition, it has to be pointed out that the IA scheme provides a simple 
physical intepretation of Nachtman's variable, whose definition was originally obtained in a 
totally different fashion 0. 

As F\ and F 2 are known to scale in x and the Bjorken variable is nearly proportional to 
y, they obviously scale in y as well. Figs. [7] and [| show that F\ and F 2 do indeed exhibit 
^-scaling at fixed |q|, whereas W 2 does not. Note that the data plotted in Figs. [7| and |8| 
span a wide range of Q 2 , extending from 1.7 to 37.5 GeV 2 . 

It is very important to realize that, while y essentially coincides with —Mx, the scaling 
function F(y) cannot be identified with either Fi or F 2 . The structure of F(y) is entirely 
dictated by the internal dynamics of the target, whereas both F 1 = Ma\F and F 2 = ua 2 F 
contain part of the cross section describing the elecron-constituent interaction. The shape of 
Fi and F 2 is in fact strongly affected by the presence of the electromagnetic cross section. For 
example, at y = gauge invariance requires F 2 to vanish, whereas F\ must be proportional 
to the photoabsorption cross section. The ^-scaling behavior displayed by Fi and F 2 is a 
consequence of the fact that, besides the proton response F, ai and va 2 also scale. This 
feature is illustrated in Fig.[| showing that ~5\ and vo 2 , plotted as a function of y, do not 
depend on |q| for |q| > 10 GeV. 

The different pictures emerging from x- and ^-scaling analyses can only be reconciled 
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making the standard assumption of parton model that the constituent mass be negligibly 
small. As shown in fig. [It], as m — > 

O 2 

W 1 ^l , ^-^p (48) 

implying in turn that, in the Bjorken limit, 

F 1 = MF(y) , F 2 = v^- F(y) = 2xF 1 , (49) 

and the standard picture of DIS is recovered. However, it has to be emphasized that, 
although in textbook derivations the requirement m ~ is often introduced as a necessary 
condition for scaling in DIS, the present analysis shows that scaling occurs irrespective of 
the constituent mass. 



V. IMPLICATIONS OF THE PROPOSED ANALYSIS 

The differences between ^/-scaling analysis and the standard x-scaling analysis, based on 
the parton model, have several interesting implications for both the theoretical description 
of DIS and the interpretation of experimental data. 

The fact that the data appear to be compatible with the naive description of the proton 
in terms of massive bound constituents suggests that past attempts to describe DIS within 
the constituent quark model (CQM) [pO| , pl"H may have to be reconsidered. 

Within the CQM, the valence quark distributions at low resolution scale Ql ~ //q ~ 
(M/3) 2 , u v (x,Ql) and d v (x,Ql), are obtained from quark momentum distributions and 
used as a starting point for QCD evolution to the Q 2 scale relevant to DIS data. As an 



illustrative example of this procedure, fig. 11 shows the valence quarks contribution to 



the proton structure function F 2 , evaluated in ref. []2tJ using a harmonic oscillator quark 
momentum distribution, before (dashed line) and after (dot-dash line) QCD evolution until 
Q 2 = 15 GeV 2 . The authors of ref. [f2(J chose a constituent mass m = M/3 and adjusted 
the harmonic oscillator frequency in such a way as to reproduce the experimental value of 
the proton rms charge radius. 
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The approach proposed in this paper provides a consistent alternative framework to 
calculate the proton structure functions from CQM momentum distributions. For any given 
momentum distribution n(p), F2 = vWi can in fact be obtained from Eq.(|3T|) using P(p) = 
n(p)5(po — Pq), p being defined in terms of the parameter B through Eqs. P5|) and (^6|). 

The structure function F2 obtained using the same constituent mass and momentum 
distribution as in ref. [2(J and B = 200 MeV is shown by the solid line in fig. [El. Note that 



in the CQM calculation of the distributions u v (x, Qq) and d v (x, Qq), the constituent quarks 
are assumed to be on mass shell. As a consequence, M-r, = M — m, implying in turn Bq = 0. 

Inclusion of binding, described by Bo 7^ 0, leads to the appearance of strength in the 
region y > — .10 GeV, where the structure function of ref. |2(J vanishes. The mechanism 



responsible for this feature can be readily understood with the help of fig. 12, showing the 



domains of the [M\, |p|) plane relevant to the calculation of the proton structure functions 



at |q| = 10 GeV and y = 0, -200 MeV and -350 MeV using Eq-flp). 

In the CQM, as the quark energy distribution is a 5-function requiring B = 0, the 
integration is carried out along the horizontal line M\ = (M — m) 2 , corresponding to 
the lower limit of the kinematically allowed domains. The maximum of the response is 
located to y — —m, because in this case the lower limit of the momentum integration is 
\Pmin\ = 0. Larger y corresponds to larger |p m j n |, leading in turn to a smaller response, as 
the quark momentum distribution is a rapidly decreasing function of |p|. For example, at 
y = —.10 GeV ^(|pmm|)/^(|p| = 0) ~ 1%. Setting B ^ 0, which amounts to carrying out 
the integration along the horizontal line M\ = (M — m + Bq) 2 , leads to a decrease of |p m j n | 
at y > —m, i.e. to an increase of the response. At y = — .1 GeV, one finds that B = and 
200 MeV correspond to |p m m| = .20 and .01 GeV, respectively. 

The fact that scaling in y occurs irrespective of the values of m and B , the resulting 
scaling function being roughly determined by the difference m — B , i.e. by the mass of the 
spectator system M-ji, poses a serious problem for the extraction of F(y) from the data and 
its interpretation. Even under the standard parton model assumption of negligibly small 
mass, allowing the constituent to be off mass shell leads to a shift of F(y) towards larger y. 
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In ref. it has been pointed out that if binding effects are large, i.e. if the energy 
distribution of the constituents extends into the region of large M^, which dominates the 
response at y ~ (see fig. [L2|), they may result in part of the response being pushed into 
the timelike region v > |q|, inaccessible to electron scattering. The occurrence of strength 
located in the timelike region is a well known feature of the response of interacting many- 
body systems, simply reflecting the fact that the there is no a priori reason preventing the 
target from being left in a final state with energy E n > Eo+\q\, where Eq is the energy of the 
initial state. On the other hand, in the case of noninteracting systems the inequality v < |q| 
is always satisfied and the response can be shown to vanish identically for y > 0. This issue 



has been recently investigated in ref. [22], whose authors have shown that ~ 10 % of the 
strength associated with the response of a confined relativistic particle is indeed located in 
the region of positive y. 

Obviously, the presence of strength in the timelike region would make the experimental 
verification of sume rules involving the proton response, e.g. the Gross-Llewellyn Smith sum 
rule [2^], impossible. This problem has long been recognized in the analysis of electron- 
nucleus scattering data. For example, integration of the measured charge response of iron 
at momentum transfer |q| = 570 MeV yields only ~ 90% of the nuclear charge Z [TM . 



VI. CONCLUSIONS 

The results described in this paper show that the scaling analysis of scattering pro- 
cesses off many-body systems, successfully employed to describe neutron-liquid helium and 
electron-nucleus scattering, can be extended to the case of DIS of electrons by protons. 

The data exhibit a remarkable scaling behavior in the natural variable of many-body 
theory, whose definition and interpretation naturally emerge from the analysis of the scat- 
tering process in the target rest frame. The onset of y-scaling follows from the assumptions 
underlying the IA picture. No further assumptions on the mass and binding energy of proton 
constituents are required. 
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While in the Bjorken limit y essentially concides with x, up to a constant factor, the 
scaling function F(y) cannot be trivially related to the structure functions Fi(x) and F 2 (x), 
whose behavior is dictated by both the internal proton dynamics and the electromagnetic 
electron-constituent interaction. The correspondence between F(y) and F\(x) and F2(x) 
can only be established in the limit of vanishing constituent mass. 

The function F(y), which can be identified with the target response in the scaling limit, 
turns out to be significantly affected by binding effects, the position of its peak being mainly 
dictated by the mass of the spectator system M-r,. 

The structure function F2 calculated from Eq . (|3l|) appreciably differs from the prediction 
of the CQM at low resolution scale. Comparison between the solid and dashed curves of fig. 



12] shows that the treatment of binding and mass according to the approach discussed in this 
paper produces an effect qualitatively similar to that of QCD evolution, shifting strength 
towards larger y. The results displayed in fig. [12] suggest that the role of QCD evolution 
should be reassessed, using the structure function represented by the solid line as a starting 
point. 

Binding effects, resulting in a shift of F(y) towards larger values of y, can also push a 
fraction of the strength into the region of positive y, inaccessible to DIS experiments. The 
proposed analysys, being carried out at fixed |q|, appears to be best suited to address this 
feature, as moving along lines of constant |q| one intersects the photon line v = |q| and 
enters the timelike region, while moving along lines of constant Q 2 > 0, as generally done 
in the analysis of DIS, the photon line is only asymptotically approached |J. 

Finally, it has to be pointed out that the observation of scaling does not necessarily 
imply the validity of the assumptions underlying the IA picture, i.e. dominance of incoher- 
ent scattering and absence of FSI. It has been shown that, both in neutron-liquid helium 



scattering [2S| and in electron-nucleus scattering ||26|| , as long as the cross section describing 
the interactions between target constituents does not depend upon |q|, FSI do not pre- 
vent the occurrence of scaling. However, in presence of FSI the scaling function extracted 
from the data can not be related to the initial momentum and energy distribution of target 
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constituents. 

The results of ref. |Z7| indicate that sizable FSI effects in the response of a scalar rel- 
ativistic particle confined by a linear potential persist in the kinematical regime of large 
|q|, in which y-scaling is clearly observed. The role of FSI in DIS and their impact on the 
interpretation of the measured proton strucure functions have been recently discussed in 



refs. [p8l,|29|l ) whose authors conclude that the DIS cross section is indeed affected by final 
state rescatterings, and can not be simply related to the parton distributions in the initial 
state. 

ACKNOWLEDGMENTS 

The author is deeply indebted to Vijay R. Pandharipande and Ingo Sick for many illu- 
minating discussions on different issues related to the subject of this paper. The hospitality 
of the Theory Group at the Thomas Jefferson National Accelerator Facility and the par- 
tial support of the U.S. Department of Energy during the completion of this work are also 
gratefully acknowledged. 

APPENDIX A: TARGET ELECTROMAGNETIC TENSOR WITHIN IA 

The general expression of the target electromagnetic tensor 

W^ = J2 I d 3 p n (0\J^n)(n\r\0)8^(P + q-P n ) (Al) 

n J 

drastically simplifies under the assumptions underlying the IA scheme, which lead to rewrite 
the target current and final state as (the index % labels the struck constituent) 

^-EiT (A2) 

i 

and 

|n) - \p',K) = \i(p , ))®\K(p R )) . (A3) 
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In the above equations jf is the current associated to a target constituent, whereas the states 
|«(p')) an d \7Z(p R )) describe the struck constituent and the spectator system, respectively 
(see the discussion following eq.(|9|)). The structure of the IA final state can be exploited to 
rewrite the sum appearing in Eq.( |Al| ) replacing 

£ \n){n\ -> / dV |p'Xp'| £ I^X^I • (A4) 

Substitution of Eqs. (|A2|) -( |A"4]) into Eq. (|Al|) and insertion of complete sets of free particle 
states, satisfying 

Jd 3 p\p)(p\ = I , (A5) 

yields the expression 

w = EE J d* P d 3 P ' \{o\ P ,K)\ 2 (pbri P ')(p'bTlp) 

x5^(q + p-p')5(v + E -E n -E pl ) , (A6) 



i n 



that can be further rewritten in a more compact form in terms of the constituent momentum 
end energy distribution P(p), defined by Eq.flT2"). As a result we obtain 



w^ = y: J d 4 P p(p) J dV(pbriP / )(p / bTlp)5 (3) (q+p-p , )^+po-^ P . (A7) 

i 

Finally, we can substitute in the above equation the definition of the electromagnetic tensor 
corresponding to scattering at four-momentum tranfer q = {y, q) off a free constituent 
carrying momentum p: 

< = J ^'(pbTipOip'bTi^^Hq+p-pO^ + ^p-^pO • (as) 

Comparison of Eqs. (|A7|) and ( |A8|) shows that W^ u can be rewritten in the form 

w ^=y: Jd 4 P p(p)wr(p,q) , (A9) 

i 

where p = (E p , p) and q = (u, q), with v = E p+({ — E p = v + p — E p . 
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APPENDIX B: RESTORATION OF GAUGE INVARIANCE 



The definition of the tensor W^ v given by Eq. obviously leads to breaking gauge 
invariance, as 



^r(M)^o (bi) 



by construction. 



The prescription proposed in ref. [12] to circumvent this problem amounts to using 



IA only for the transverse and time components of the current jf. Denoting by Wa V the 



electromagnetic tensor of ref. |12| and choosing the z-axis along the direction of q we can 
then write 

wT(p, q) = Wi v (p, o) , t*, v ^ 3 ■ ( B2 ) 

The longitudinal current jf is defined in terms of the time component j® in such a way as 
to satisfy the continuity equation. As a result, the corresponding elements of the tensor w^ v 
read 

wf{p, q)= f y^j J wf(p, q) (B3) 



and 



wf(p,q)= (j^j w°°(p,q) . (B4) 



The above procedure to reconcile the IA scheme with gauge invariance is manifestly 
non unique. However, it has to be emphasized that, since according to ref. fl2| the cur- 
rent satisfies the relation j 3 = (z//|q|)j°, while the IA current satisfies j 3 = (y/\q\)j°, the 
difference 

~f-f = ^ = ^l, (B5) 

|q| |q| 

providing a measure of the violatation of gauge invariance, becomes vanshingly small in the 
|q| — > oo limit. 
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FIG. 1. Scaling functions F(y), denned as in Eq, (|l9|) , measured by neutron scattering off 
superfluid 4 He at T=1.6 °K |]. The data sets are labelled according to the magnitude of the 
momentum transfer |q|. 
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FIG. 2. Cross sections for scattering of 4 GeV electrons off Iron, plotted as a function of the 
electron energy loss. The data sets are labelled according to the momentum transfer at the quas 
ielastic peak [fbl. 
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FIG. 3. Scaling functions obtained from the definitions of Eqs.(|35|) and (|36|) using the cross 



sections shown in fig.g |14]. 
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FIG. 11. Valence quarks contribution to the proton F2, plotted as a function of y at constant 
|q| = 10 GeV. The dashed and dash-dot lines show the results of the constituent quark model of 
ref. [20 1 before and after QCD evolution until Q 2 = 15 GeV 2 , respectively. The solid line shows F2 
obtained within the approach described in this paper, using the same quark mass and momentum 



distribution as in ref. [§0| and B = 200 MeV 
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FIG. 12. Domains of the (M^, |p|) plane relevant to the calculation of the proton structure 
functions at |q| = 10 GeV and y = (solid line), —200 MeV (dashed line) and —350 MeV (dot-dash 
line) from Eq.(pTI). The results have been obtained using m = M/3. 



34 



